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We investigate the correlation structure of pure A''-mode Gaussian resources which can be experimentally 
generated by means of squeezers and beam splitters, whose entanglement properties are generic. We show that 
those states are specified (up to local unitaries) by N{N — l)/2 parameters, corresponding to the two-point 
correlations between any pair of modes. Our construction yields a practical scheme to engineer such generic- 
entangled A'^-mode pure Gaussian states by linear optics. We discuss our findings in the framework of Gaussian 
matrix product states of harmonic lattices, raising connections with entanglement frustration and the entropic 
area law. 
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Introduction. — Multipartite entanglement in pure states 
of many systems is a founding property and a crucial resource 
for quantum information science, yet its complete theoretical 
understanding is still lacking. A basic property of entangle- 
ment is its invariance under unitary operations performed lo- 
cally on the subsystems. To describe entanglement efficiently, 
is thus natural to lighten quantum systems of the unneces- 
sary degrees of freedom adjustable by local unitaries (LUs), 
and to classify states according to standard forms representa- 
tive of LU equivalence classes 1 1|. Alongside the traditional 
qubit-based approach, quantum information with continuous 
variables (CV) is a burgeoning field mainly spinning around 
the theory and applications of entanglement in Gaussian states 

i. 

In this Letter we address the question of how many phys- 
ical resources are really needed to engineer and characterize 
entanglement in pure Gaussian states of an arbitrary number 
of modes, up to LU operations. For states of < 3 modes, 
it has been shown that such a number of minimal degrees of 
freedom scales as N{N — l)/2 QQ]. For a higher number 
of modes, however, a richer structure is achievable by pure 
Gaussian states, as from the normal form of Ref. |5] a min- 
imal number of parameters given by N{N — 2) can be in- 
ferred |6]. A random state of A^ > 4 modes, selected ac- 
cording to the uniform distribution over pure Gaussian states, 
will be thus reducible to a form characterized by such a num- 
ber of independent quantities. However, in practical realiza- 
tions of CV quantum information one is interested in states 
which, once prepared with efficient resources, still achieve an 
almost complete structural variety in their multipartite entan- 
glement properties. Such states will be said to possess generic 
entanglement 0], where generic means practically equivalent 
to that of random states, but engineered (and described) with 
a considerably smaller number of degrees of freedom. 

Precisely, we define as "generic-entangled" those Gaussian 
states whose local entropies of entanglement in any single 
mode are independent, and bipartite entanglements between 
any pair of modes are unconstrained. Having a standard form 
for such A^-mode Gaussian states, may be in fact extremely 



helpful in understanding and quantifying multipartite CV en- 
tanglement, in particular from the theoretical point of view of 
entanglement sharing and monogamy constraints and 
from a more pragmatical approach centered on using entan- 
glement as a resource. We show that, to achieve generic en- 
tanglement, for the global pure A^-mode Gaussian state it is 
enough to be described by a minimal number of parameters 
(corresponding to the LU invariant degrees of freedom) equal 
to A^(A^ — l)/2 for any A^, and thus much smaller than the 
2N{2N + l)/2 of a completely general covariance matrix. 
Therefore, generic entanglement appears in states which are 
highly not 'generic' in the sense usually attributed to the term, 
i.e. randomly picked. Crucially, we demonstrate that generic- 
entangled Gaussian states coincide with the resources typi- 
cally employed in experimental realizations of CV quantum 
information |2], and we provide an optimal scheme for their 
state engineering. 

Preliminaries. — We consider a CV system consisting of 
A^ canonical bosonic modes, and described by the vector X ~ 
{qi,Pi, q2,P2, ■ ■ ■ An ,Pn} the field quadrature operators, 
which satisfy the commutation relations [X^, A"j] = 2iily, 
with the symplectic form O — and u) = ( q)- Gaus- 
sian states (such as vacua, coherent, and squeezed states) are 
defined by having a Gaussian characteristic function in phase 
space 101 ■ They are fully characterized by the first statistical 
moments (arbitrarily adjustable by LUs: we will set them to 
zero) and by the 2A^ x 2A^ covariance matrix (CM) cr of the 
second moments (Ty ^ ({X^, Aj})/2 - {X.i){Xj). The CM 
cr of an arbitrary A^-mode Gaussian state can be written as 
follows in terms of 2 x 2 submatrices 




Symplectic operations {i.e. belonging to the group 
Spi2N.M) = {5' e SL{2N,Wi) : S^nS = n}) act- 
ing by congruence on CMs in phase space, amount to unitary 
operations on density matrices in Hilbert space. Any A^-mode 
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Gaussian state can be transformed by symplectic operations 
in its Williamson diagonal form ||9J i^, such that cr = S^uS, 
with 1/ — diagjz^i, 1^1, . . . i^TV, j/at}. The quantities > 1 
are the symplectic eigenvalues of cr 1 10|. 

We define the symplectic rank H of a CM as the number 
of its symplectic eigenvalues different from 1, corresponding 
to the number of non-vacua normal modes. Any pure state 
has N = 0. The CM ctP of any A^-mode pure Gaussian state 
satisfies the matrix identity 01 — n <tP r2 cr^ = 1. As a con- 
sequence of the Schmidt decomposition, applied at the CM 
level 1 1 1 J for the bipartition — + N), any 

{N — l)-mode reduced CM of the CM cr^ of a A^-mode pure 
Gaussian state has symplectic rank N = 1. 

Minimal number of parameters. — Adopting the above 
definition of generic entanglement, we prove now the main 
Proposition 1: A generic-entangled A^-mode pure Gaussian 
state is described, up to local symplectic (unitary) operations, 
by N{N — l)/2 independent parameters. 
Proof. Let us start with a A^-mode pure state, described by a 
CM a-P = cr with all single-mode blocks in diagonal form: 
we can always achieve this by local single-mode Williamson 
diagonalizations in each of the A^ modes. Let cr^^ = (T2^... jv 
be the reduced CM of modes (2, . . . , A^). It can be diagonal- 
ized by means of a symplectic S2,--- ,n, and brought thus to its 
Williamson normal form, characterized by a symplectic spec- 
trum {a,l,-- - ,1}, where a = -/Detal. Transforming cr by 
S = S2,--- ,N, brings the CM into its Schmidt form, con- 
stituted by a two-mode squeezed state between modes 1 and 2 
(with squeezing a), plus N — 2 vacua 1 11, 12il . 

All A^-mode pure Gaussian states are thus completely spec- 
ified by the symplectic S'2,... .n, plus the parameter a. Alter- 
natively, the numbers of parameters of cr is also equal to those 
characterizing an arbitrary mixed A^ — 1 Gaussian CM, with 
symplectic rank H = 1 (i.e. with N — 2 symplectic eigenval- 
ues equal to 1). This means that, assigning the reduced state 
cr2,- ,N, we have provided a complete description of cr. In 
fact, the parameter a is determined as the square root of the 
determimant of the CM 0-2 ... jy. 

We are now left to compute the minimal set of parameters of 
an arbitrary mixed state of A^ — 1 modes, with symplectic rank 
H = 1. While we know that for A^ > 4 this number is equal to 
N{N — 2) in general j^, we want to prove that for generic- 
entangled Gaussian resource states this number reduces to 

En = N{N -l)/2. (2) 

We prove it by induction. For a pure state of one mode only, 
there are no reduced "zero-mode" states, so the number is 
zero. For a pure state of two modes, an arbitrary one-mode 
mixed CM with H = 1 is completely determined by its own 
determinant, so the number is one. This shows that our law 
for Sjv holds true for A^ = 1 and N ^ 2. 

Let us now suppose that it holds for a generic A^, i.e. we 
have that a mixed (A^ — l)-mode CM with N = 1 can be 
put in a standard form specified by A^(A^ — l)/2 parameters. 
Now let us check what happens for a (A^ + l)-mode pure state, 
i.e. for the reduced A^-mode mixed state with symplectic rank 



equal to 1. A general way (up to LUs) of constructing a A'^- 
mode CM with H = 1 yielding generic entanglement is the 
following: (a) take a generic-entangled (A^ — l)-mode CM 
with H = 1 in standard form; (b) append an ancillary mode 
(cr^v) initially in the vacuum state (the mode cannot be ther- 
mal as H must be preserved); (c) squeeze mode A^ with an 
arbitrary s (one has this freedom because it is a local sym- 
plectic operation); (d) let mode A^ interact couplewise with all 
the other modes, via a chain of beam-splitters fTill with arbi- 
trary transmittivities bi,N, with i = 1, • • • , A^ — 1 lll4ll : (e) if 
desired, terminate with A^ suitable single-mode squeezing op- 
erations (but with all squeezings now fixed by the respective 
reduced CM's elements) to symplectically diagonalize each 
single-mode CM. 

With these steps one is able to construct a mixed state 
of A^ modes, with the desired rank, and with generic (LU- 
invariant) properties for each single-mode individual CM. We 
will show in the following that in the considered states the 
pairwise quantum correlations between any two modes are un- 
constrained. To conclude, let us observe that the constructed 
generic-entangled state is specified by a number of parame- 
ters equal to: N{N — l)/2 (the parameters of the starting 
(A^— l)-mode mixed state of the same form) plus 1 (the initial 
squeezing of mode A^) plus A^ — 1 (the two-mode beamsplitter 
interactions between mode A^ and each of the others). Total: 
iN + l)N/2 = EN+i. ■ 

Quantum state engineering. — Following the ideas of 
the above proof, a physically insightful scheme to produce 
generic-entangled A^-mode pure Gaussian states can be read- 
ily presented (see Fig. 0. It consists of basically two main 
steps: (1) creation of the state in the 1|(A^ — 1) Schmidt de- 
composition; (2) addition of modes and entangling operations 
fl3ll between them. One starts with a chain of A^ vacua. 

First of all (step 1), the recipe is to squeeze mode 1 of an 
amount s, and mode 2 of an amount 1/s (i.e. one squeezes 
the first mode in one quadrature and the second, of the same 
amount, in the orthogonal quadrature); then one lets the two 
modes interfere at a 50 : 50 beam splitter One has so created 
a two-mode squeezed state between modes 1 and 2, which 
corresponds to the Schmidt form of cr with respect to the 
1|(A^ — 1) bipartition. The second step basically corresponds 
to create the most general mixed state with H = 1, of modes 
2, • • • , A^, out of its Williamson diagonal form. This task can 
be obtained, as already sketched in the above proof, by letting 
each additional mode interact step-by-step with all the previ- 
ous ones. Starting with mode 3 (which was in the vacuum like 
all the subsequent ones), one thus squeezes it (of an amount 
ra) and combines it with mode 2 via a beam-splitter (charac- 
terized by a transmittivity 62,3)- Then one squeezes mode 4 
by r4 and lets it interfere sequentially both with mode 2 (with 
transmittivity 62,4) and with mode 3 (with transmittivity 63,4). 
This process can be iterated for each other mode, as shown 
in Fig. n until the last mode A^ is squeezed (r^) and entan- 
gled with the previous ones via beam-splitters with respective 
transmittivities bi^N, i — 2, ■ ■ ■ , N — 1. Step 2 describes the 
distribution of the two-mode entanglement created in step 1, 
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FIG. 1: (Color online) How to create a generic-entangled A'^-mode 
pure Gaussian state. White balls are vacua, while each color de- 
picts a different single-mode determinant {i.e. different degrees of 
local mixedness). Vertical arrows denote single-mode squeezing op- 
erations, while horizontal circle-ended lines denote beam-splitting 
operations j between modes i and j. See text for details. 



among all modes. 

The presented prescription enables to create a generic form 
(up to LUs) of multipartite entanglement among N modes 
in a pure Gaussian state, by means of active (squeezers) and 
passive (beam-splitters) linear optical elements. What is rel- 
evant for practical applications, is that the state engineer- 
ing is implemented with minimal resources. Namely, the 



process is characterized by one squeezing degree (step 1), 
plus N — 2 individual squeezings for step 2, together with 
i = {N - 1){N - 2)/2 beam-splitter transmittivities, 
which amount to a total of N{N — l)/2 = S^v quantities. 
The optimally produced Gaussian states can be readil y im ple- 
mented for iV-party CV communication networks O llST . 

Standard form. — The special subset of pure A^-mode 
Gaussian states emerging from our constructive proof exhibits 
a distinct property: all correlations between "position" qi and 
"momentum" pj operators are vanishing. Looking at Eq. Q, 
this means that such a generic-entangled pure Gaussian state 
can be put in a standard form where all the 2x2 subma- 
trices of its CM are diagonal. The diagonal subblocks cTi 
can be additionally made proportional to the identity by lo- 
cal Williamson diagonalizations in the individual modes. This 
standard form for generic-entangled TV-mode Gaussian states, 
as already mentioned, can be achieved by all pure Gaussian 
states for /V 2 |31 and /V = 3 |4]; for > 4, pure Gaus- 
sian states can exist whose number of independent parameters 
scales as N{N — 2) f3\ and which cannot thus be brought 
in the q-p block-diagonal form. Interestingly, all pure Gaus- 
sian states in our considered block-diagonal standard form are 
ground states of quadratic Hamiltonians with spring-like in- 
teractions 1 16|. 

Vanishing q-p covariances imply that the CM can be writ- 
ten as a direct sum cr'P — Vq ® Vp, when the canonical 
operators are arranged as {qi, . . . , qjv, Pi, ■ ■ ■ ,Pn}- More- 
over, the global purity of ctP imposes Vp = Vq^. Named 
(yQ)ij = '"Qij and (yp)hk = wPfefc, this means that 
each vpf^^ is a function of the {wq^ }'s. The additional N 
Williamson conditions vp.. = vq.. fix the diagonal elements 
of Vq. The standard form is thus completely specified by the 
off-diagonal elements of the symmetric N x N matrix Vq, 
which are, as expected, N{N — l)/2 = S^r. Proposition 
1 acquires now a remarkable physical insight: the structural 
properties of the generic-entangled A^-mode Gaussian states, 
and in particular their bipartite and multipartite entanglement, 
are completely specified (up to LUs) by the 'two-point cor- 
relations' VQ.. = (qiqj) between any pair of modes. For in- 
stance, the entropy of entanglement between one mode (say 
i) and the remaining — 1 modes, which is monotonic in 
Det CTi 1 10], is completely specified by assigning all the pair- 
wise correlations between mode i and any other mode j ^ i, 
as Dct (T j = 1 — '^j-ii Dct Sij . The rationale is that entangle- 
ment in such states is basically reducible to a mode-to-mode 
one. This statement, strictly speaking true only for the pure 
Gaussian states for which Proposition 1 holds, acquires a gen- 
eral validity in the context of the modewise decomposition of 
arbitrary pure Gaussian states liTllfl2ll . This correlation pic- 
ture breaks down for mixed Gaussian states, where also clas- 
sical, statistical-like correlations arise. 

Gaussian matrix product states. — As an application, let 
us consider Gaussian matrix product states (GMPS), defined 
as A^-mode states obtained by taking a fixed number, M, of 
infinitely entangled ancillary bonds [Einstein-Podolski-Rosen 
(EPR) pairs] shared by adjacent sites, and applying an arbi- 
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trary 2M 1 Gaussian operation "Pl'l on each site i = 
1, . . . , N. The projections T'I'l can be described in terms of 
isomorphic (2M + l)-mode pure Gaussian states with CM 
7W, the building blocks frT]. 

It is conjectured that all pure iV-mode Gaussian states can 
be described as GMPS. Here we provide a lower bound on the 
number AI of ancillary bonds required to accomplish this task, 
as a function of N. We restrict to ground states of harmonic 
chains with spring-like interactions. With a simple counting 
argument, the total number of parameters of the initial chain 
of building blocks should be at least equal to that of the target 
state, i.e. N{2M + l)(2A f)/2 > N{N - l)/2 which means 
M > IntPart[(V4^ - 3 - l)/4]. This impUes, for instance, 
that to describe generic states with at least N > 7 modes, 
a single EPR bond per site is no more enough (even though 
the simplest case of M = 1 yields interesting families of iV- 
mode GMPS for any TV LJ). The minimum M scales as 
N^^^, diverging in the field limit ^ 00. As infinitely many 
bonds would be necessary (and maybe not even sufficient) to 
describe generic infinite harmonic chains, the matrix product 
formalism is probably not helpful to prove or disprove area 
law statements for critical bosons (complementing the known 
results for the non-critical case 1 18]), which in general do not 
fall in special subclasses of finite-bonded GMPS. 

The matrix product picture however effectively captures 
the entanglement distribution in translationally invariant iV- 
mode harmonic rings JItIi . In this case the GMPS build- 
ing blocks are equal at all sites, -y^'^^ = 7 Vz, while the 
number of parameters Eq. (|2jl of the target state reduces to 
the number of independent pairwise correlations (only func- 
tions of the distance between the two sites), which by our 
counting argument is Oat = {N — N mod 2)/2. The cor- 
responding threshold for a GMPS representation becomes 
M > IntPart[(V8e^vTT - l)/4]. As 9^ is bigger for 
even N, so it is the resulting threshold, which means that in 
general a higher number of EPR bonds is needed, and so more 
entanglement is inputed in the GMPS projectors and gets dis- 
tributed in the target iV-mode Gaussian state, as opposed to 
the case of an odd N. This clarifies why nearest-neighbour 
entanglement in ground states of pure translationally invariant 
A^-mode harmonic rings (which belong to the class of states 
characterized by Proposition 1) is frustrated for odd N Q^. 

Conclusions. — In this Letter we studied pure A^-mode 
Gaussian states of CV systems, aimed to eradicate the mini- 
mal number of degrees of freedom responsible for the generic 



nonlocal features of the states. We showed that a crucial 
subclass of such states, employed as typical resources in CV 
quantum information, can be put in a standard form described 
(up to local unitaries) by N{N — l)/2 parameters, corre- 
sponding to the ground states of a quadratic Hamiltonian with 
spring-like couplings. This form encompasses all pure Gaus- 
sian states for < 3. We operationally related these pa- 
rameters with the active and passive transformations needed 
to prepare the state. In general we interpreted those degrees 
of freedom as the two-point correlations between any pair of 
modes, which are thus responsible for the structure of generic 
entanglement in Gaussian states. It would be worth to investi- 
gate the exotic properties of multipartite entanglement arising 
in Gaussian states which cannot be brought in the standard 
form described here for A^ > 4. 
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